Dynamic Resonance of Light in Fabry-Perot Cavities 
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The dynamics of light in Fabry-Perot cavities with varying length and input laser frequency are 
analyzed and the exact condition for resonance is derived. This dynamic resonance depends on the 
light transit time in the cavity and the Doppler effect due to the mirror motions. The response of 
the cavity to length variations is very different from its response to laser frequency variations. If the 
frequency of these variations is equal to multiples of the cavity free spectral range, the response to 
length is maximized while the response to the laser frequency is zero. Implications of these results 
for the detection of gravitational waves using kilometer-scale Fabry-Perot cavities are discussed. 
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Fabry-Perot cavities, optical resonators, are commonly 
utilized for high-precision frequency and distance mea- 
surements 0. Currently, kilometer-scale Fabry-Perot 
cavities with suspended mirrors are being employed in 
efforts to detect cosmic gravitational waves ||, ^. This 
apphcation has stimulated renewed interest in cavities 
with moving mirrors ^ 01 and motivated efforts 
to model the dynamics of such cavities on the computer 
^, |l^, |ri|, |l^ . Recently, several studies addressed the 



process of lock acquisition in which the cavity mirrors 
move through the resonance positions |l3|, Q . In this 
process, the Doppler effect due to the mirror motions 
impedes constructive interference of light in the cavity 
giving rise to complex field dynamics. In contrast, Fabry- 
Perot cavities held in the state of resonance have usually 
been treated as essentially static. In this letter, we show 
that resonant cavities also have complex field dynamics 
and we derive the condition for dynamic resonance. Our 
analysis is developed for the very long Fabry-Perot cav- 
ities of gravitational wave detectors, but the results are 
general and apply to any cavities, especially when the fre- 
quencies of interest are close to the cavity free spectral 
range. 

We consider a Fabry-Perot cavity with a laser field in- 
cident from one side as shown in Fig. |l|. Variations in the 
cavity length are due to the mirror displacements Xa{t) 
and Xb{t) which are measured with respect to reference 
planes a and h. The nominal light transit time and the 
free spectral range (FSR) of the cavity are defined by 



FIG. 1: Mirror positions and fields in a Fabry-Perot cavity. 



fine the slowly-varying field as E{t) — £{t) exp{—iujt). 
To properly account for the phases of the propagating 
fields, their complex amplitudes are defined at fixed lo- 
cations, reference planes ai and 02, as shown in Fig. |^. 
(The small offset e is introduced for convenience and can 
be set to zero at the end of calculations.) 

The equations for fields in the cavity can be obtained 
by tracing a wavefront during its complete round-trip in 
the cavity (starting from the reference plane 02). The 
propagation delays ti and T2 depend on the mirror posi- 
tions and are given by 



CTi = L — e + xi,{t — Ti), 

CT2 = e - Xa{t - T2). 



Then the fields in the cavity satisfy the equations: 



T = L/c, 



11 /T. 



(1) 



E'{t) = -rb£;(t-2Ti)e-2--i^ 
E{t) = ^raE'{t-2T2)e 



— 2iuJT2 



+ 



The field incident upon the cavity and the field circu- 
lating in the cavity are described by plane waves with 
nominal frequency w and wavenumber k (fc = uj/c). Vari- 
ations in the laser frequency are denoted by Suj{t). We 
assume that the mirror displacements are much less than 
the nominal cavity length and that the deviations of the 
laser frequency are much less than the nominal frequency. 

At any given place the electric field £ in the cavity os- 
cillates at a very high frequency: £{t) oc exp(iujt). For 
simplicity, we suppress the fast-oscillating factor and de- 



taEinit-2e/c), 



(2) 
(3) 



(4) 
(5) 



where and ri, are the mirror reflectivities, and ta is the 
transmissivity of the front mirror. 

Because the field amplitudes E and E' do not change 
significantly over times of order Xa,b/c, the small varia- 
tions in these amplitudes during the changes in propaga- 
tion times due to mirror displacements can be neglected. 
Furthermore, the reference planes a and b can be chosen 
so that the nominal length of the Fabry-Perot cavity be- 
comes an integer multiple of the laser wavelength, making 



2 



exp{—2ikL) = 1. Finally, the offset e can be set to zero, 
and Eqs. (^-(^) can be combined yielding one equation 
for the cavity field 

E{t) = taEinit) + ranE{t - 2T) exp[~2ikSL{t)]. (6) 

Here SL{t) is the variation in the cavity length "seen" by 
the light circulating in the cavity. 



SL{t) ^Xb{t~T)-Xa{t). 



(7) 



Note that the time delay appears in the coordinate of 
the end mirror, but not the front mirror. This is simply 
a consequence of our placement of the laser source; the 
light that enters the cavity reflects from the end mirror 
first and then the front mirror. For 6L = 0, Laplace 
transformation of both sides of Eq. (^) yields the basic 
cavity response function 



H{s) 



tn 



Ejs) ^ 

Einis) I -Tar be 



-2sT ' 



(8) 



where tildes denote Laplace transforms. 

The static solution of Eq. (|^) is found by considering a 
cavity with fixed length (6L — const) and an input laser 
field with fixed amplitude and frequency {A, 6uj = const). 
In this case the input laser field and the cavity field are 
given by 



E{t) = Eo 



iSujt 



where Eq is the amplitude of the cavity field, 

taA 



Eo 



1 — rarb exp[— 2i(r Sw + k SL)] 



(9) 
(10) 



(11) 



The cavity field is maximized when the length and the 
laser frequency are adjusted so that 



Slu SL 

LU L 



(12) 



This is the well-known static resonance condition. The 
maximum amplitude of the cavity field is given by 



E 



taA 



1 - ran 



(13) 



Light can also resonate in a Fabry-Perot cavity when 
its length and the laser frequency are changing. For a 
fixed amplitude and variable phase, the input laser field 
can be written as 

Ein{t) ^ A e"''^*\ (14) 
where (j){t) is the phase due to frequency variations. 



Jo 



(15) 



Then the steady-state solution of Eq. is 

E{t) ^Ee"*'^'\ 



(16) 



where the amplitude E is given by Eq. (13) and the phase 
satisfies the condition 



(j){t) - (l){t - 2T) = -2k 5L{t). 



(17) 



Thus resonance occurs when the phase of the input laser 
field is corrected to compensate for the changes in the 
cavity length due to the mirror motions. The associated 
laser frequency correction is equal to the Doppler shift 
caused by reflection from the moving mirrors 

5uj(t) - 5uj{t - 2T) = -2—uj, (18) 

c 

where v{t) is the relative mirror velocity {v = dSL/dt). 
The equivalent formula in the Laplace domain is 



C(s 



SQ{s) SL{s) 



(19) 



where C(s) is the normalized frcquency-to-length trans- 
fer function which is given by 



Cis) 



1 



~2sT 



2sT 



(20) 



Eq. ( |l9|) is the condition for dynamic resonance. It must 
be satisfied in order for light to resonate in the cavity 
when the cavity length and the laser frequency are chang- 
ing. 

The transfer function C(s) has zeros at multiples of 
the cavity free spectral range. 



(21) 



where n is integer, and therefore can be written as the 
infinite product. 



C{s) 



-sT 



n 



(22) 



which is useful for control system design Q. 

To maintain resonance, changes in the cavity length 
must be compensated by changes in the laser frequency 
according to Eq. (^ . If the frequency of such changes is 
much less than the cavity free spectral range, C(s) ~ 1 
and Eq. (|l^) reduces to the quasi-static approximation. 



Soj{s) SLis) 



(23) 



At frequencies above the cavity free spectral range, 
C(s) oc 1/s and increasingly larger laser frequency 
changes are required to compensate for cavity length vari- 
ations. Moreover, at multiples of the FSR, C(s) =0 and 
no frequency-to-length compensation is possible. 




FIG. 2: Bode plot of for the LIGO 4-km Fabry- Perot 

cavities. The peaks occur at multiples of the FSR (37.5 kHz) 
and their half-width (91 Hz) is equal to the inverse of the 
cavity storage time. 



In practice, Fabry-Perot cavities always deviate from 
resonance, and a negative- feedback control system is em- 
ployed to reduce the deviations. For small deviations 
from resonance, the cavity field can be described as 



E{t) = [E + dE{t)]e 



i4>(t) 



(24) 



where E is the maximum field given by Eq.(p^), and 
5E is a small perturbation {\5E\ <C |-^|). Substituting 
this equation into Eq. (^, we see that the perturbation 
evolves in time according to 

5E{t) - rarb5E{t - 2T) = 

-iranE [(/)(t) - <j){t - 2T) + 2k SL(t)] . (25) 

This equation is easily solved in the Laplace domain, 
yielding 



SE{s) = -iranE 



1 



-2sT^ 



(s) -I- 2k 5L{s) 



1 - TaThe 



-2sT 



(26) 



Deviations of the cavity field from its maximum value 
can be measured by the Pound-Drever-Hall (PDH) er- 
ror signal which is widely utilized for feedback control of 
Fabry-Perot cavities The PDH signal is obtained 

by coherent detection of phase-modulated light reflected 
by the cavity. With the appropriate choice of the de- 
modulation phase, the PDH signal is proportional to the 
imaginary part of the cavity field (Eq. (^)) and therefore 
can be written as 



SV{s) oc H{s) 



SL{s) 



Cis) 



Soj{s) 



(27) 
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FIG. 3: Bode plot of Hui{iO.). The sharp features are due to 
the zero-pole pairs at multiples of the FSR. 



error signal toward the null point, SV = 0, thus main- 
taining dynamic resonance accor ding to Eq. (|l|). 

The response of the PDH signal to either length or laser 
frequency deviations can be found from Eq. (p^). The 
normalized length-to-signal transfer function is given by 



Hl{-s) = 



1 - ran 



His) 

H{0) 1 - rane 



-2sT ■ 



(28) 



A Bode plot (magnitude and phase) of Hl is shown in 
Fig. I for the LIGO [g Fabry-Perot cavities with L = 4 
km, Ta = 0.985, and = 1. The magnitude of the 
transfer function, 



\HLm\ 



1 



Vl + sin^ ^IT 



(29) 



is the square-root of the well-known Airy function with 
the coefficient of finesse F = Aran/{^ — Tarbf'- (In op- 
tics, the Airy function describes the intensity profile of a 
Fabry-Perot cavity |^ .) 

The transfer function has an infinite number of 
poles: 



1 



Pn 



(30) 



where n is integer, and t is the storage time of the cavity, 

2T 

r^r^, W (31) 

I \n(ran)\ 

Therefore, H]^ can be written as the infinite product. 



Pn 



Pn - S 



(32) 



where H{s) is given by Eq. (g). In the presence of length 
and frequency variations, feedback control will drive the 



which can be truncated to a finite number of terms for 
use in control system design. 
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The response of a Fabry-Perot cavity to laser frequency 
variations is very different from its response to length 
variations. The normalized frequency-to-signal transfer 
function is given by 



H^{s) = C{s)Hl{s), 



or, more explicitly as 



2sT 



1 - ran 



1 - rarbe 



-2sT 



(33) 



(34) 



A Bode plot of H^, calculated for the same parameters 



The transfer function 
with n ^ 0, and poles 



as for H]^, is shown in Fig. 

has zeros given by Eq. (Ell) 
given by Eq. (|3C|). The poles and zeros come in pairs 
except for the lowest order pole, po, which does not have 
a matching zero. Therefore, can be written as the 
infinite product. 



Po 



Po 



n 



I - s/Zn 
1 - s/pn 



(35) 



where the prime indicates that n = term is omitted 
from the product. 

The zeros in the transfer function indicate that the 
cavity does not respond {SE — 0) to laser frequency de- 
viations if these deviations occur at multiples of the cav- 
ity FSR. In this case, the amplitude of the circulating 
field is constant while the phase of the circulating field is 
changing with the phase of the input laser field. 

In summary, we have shown that resonance can be 
maintained in a Fabry-Perot cavity even when the cav- 
ity length and laser frequency are changing. In this dy- 
namic resonance state, changes in the laser frequency 
and changes in the cavity length play very different roles 
(Eq. (|l9|)) in contrast to the quasi-static resonance state 
where they appear equally (Eq. (|2^)). Maintenance of 
dynamic resonance requires that the frequency-to-lcngth 
transfer function, C(s), be taken into account when com- 
pensating for length variations by frequency changes and 
vice versa. Compensation for length variations by fre- 
quency changes becomes increasingly more difficult at 
frequencies above the FSR, and impossible at multiples 
of the FSR. 

As can be seen in Fig. ^, the response of the PDH error 
signal to laser frequency variations decreases as l/fi for 
f2 ^ r^^ and becomes strongly suppressed at frequencies 
equal to multiples of the cavity FSR. In contrast, the 
response to length variations is a periodic function of 
frequency as shown in Fig. ^j. For ^ t^^, it also 
decreases as l/fi but only to the level of (1 -I- F)~2 and 
then returns to its maximum value. Thus, at multiples of 
the FSR, the sensitivity to length variations is maximum 
while the sensitivity to frequency variations is minimum. 



Both these features suggest searches for gravitational 
waves at frequencies near multiples of the FSR. However, 
because gravitational waves interact with the light as well 
as the mirrors, the response of an optimally-oriented in- 
terferometer is equivalent to Hc^{s) and not to Hl{s) 
Thus, an optimally-oriented interferometer does not 
respond to gravitational wave at multiples of the FSR. 
However, for other orientations gravitational waves can 
be detected with enhanced sensitivity at multiples of the 
cavity FSR 0- 
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